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a ATTENUATION OF STRESS WAVES IN BI- LINEAR MATERIALS | 
By Richard Skalak, 1M E, and Paul Weidlinger,2 2 F, ASCE 


_ pressure in a bi-linear medium is investigated, to obtain an approximate 
>: solution to the problem of a medium with a stress strain diagram of posi- i. 
tive curvative. One-dimensional wave equations are derived, andit is demon- 
_ strated that, beyond a given distance from the surface, the intensity of the 
peak stress and particle velocity depends only on a single parameter. This 


S result is of significance in establishing the relevant physical characteristics © 
granular soils which are subjected to intensity surface blast 


‘pressures, 


_ Granular or porous n materials, such as sand 1 and/or certain ‘types of rock, 


behavior, | characterized by a stress-strain curve (Fig. 1), of positive curva- _ 
ture, such that aa2/ae2 > 0. On unloading, typically, a a path, such 


_ The attenuation of blast pressure is of major importance in investigations 


onuntens with the effects of nuclear surface bursts, or underground structures 


a Note.—Discussion open until November 1, 1961. To extend the closing date one month, 
-a written request must be filed with the Executive Secretary, ASCE, This paper is part | 
of the copyrighted Journal of the Engineering Mechanics Division, Proceedings of the | 
American Society of Civil Engineers, Vol, 87, No. EM3, June, 1961, 

1 assoc. Prof. of Civ. Engrg., Columbia Univ., New York, (on leave 2 at Trinity Col- 

lege, Cambridge Univ., England), 

Cons, Engr., New York, N, 


‘subjected to high intensity compression, show anon- linear, strain- hardening 
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‘ and the detection of ‘underground nuclear explosions. In order to discover the 
influence of the properties of the medium onthese phenomena, the attenuation 
- of stress waves as a functionofthe parameters which define the stress-strain 
-—sTn the present paper, one- dimensional wave propagation with finite strain 
_ will be considered in a oom infinite space, , which on its free surface is sub- 


jected t toa pressure pulse of the form shown in Fig. 2. This problem 
solved by the method of characteristics, 8, 4 4 for. arbitrary stress- -Strain curves, 


“Mathematical Theory of of Compressible ‘Fluid Flow, ” by R. Von Mises, 


Press, New York, 1958. 
“Supersonic Flow and Shock Waves,” by R. Courant and K.O, Friedrich, Interscience _ 


4 
— 

FIG, 2,—TYPICAL PRESSURE LOADING 


and num numerical solutions have been obtained for the case ase of a a constant boundary _ . 
velocity. 5 Such procedure, however, requires extensive numerical integration 
and does not readily permit the assessment of the influence of the various 
‘parameters of the problem, More direct solutions tothis problem are obtain- 
able if a simplified bi-linear model of the stress-strain curve, such as shown 
Fig. 3, On this strain curve, up to the stress o¢ an 


<E= 00 (LOCKING) 


‘FIG, 3.—BI- -LINEAR STRESS- STRAI 


elastic be behavior is and the cannot be , distinguished an 
elastic material. If the medium is compressed beyond o , the path OA'A is 


- followed, but unloading occurs on the line AB, resulting in a Y permanent strain 
€¢. The medium is described by four parameters: the critical stress 0 ¢, the 


_ 5 “Propagation of Plane Waves in Granular Mate rials, Part I,” by S Katz and T. J, 
Ahrens, Rensselaer Polytechnic Inst, No, 1, Contract DA-20-115- 509-ORD- 1009, 
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permanent strain Ries and the slope a of the second branch | of the curve, which 

_ defines the elastic modulus E = tan a, and if the initial density p of the medium - 


is also given, the material is defined. 
ae now the limiting values of the parameters, the special case 
is recognized, that is, the stress-strain diagram becomes that — 


material becomes incompressible, when the 
= 0, 


ane, 196 ~ 


when E = 
«a a locking medium, 5 


with an elastic ‘modulus E= = ‘The main purpose ose of this p: paper is to in- 
vestigate variations of the compressibility within the limits LS EX %, onal 


Notation.—The letter ‘symbols adapted for use in this paper are ‘defined | 


where” they first appear, int the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, in “Appendix II. = 


WAVE EQUATIONS 


Consider fi first the case in which Oc in Fig. 3 is zero so that the stress- 
strain diagram consists of lines OB and BA. (Results for the case Oc > 0 are 


easily obtained, as will be shown later, once the case ¢ = 0 is solved. 
‘The stress wave due toa loading of the form of Fig. 2 consists of a sharp | 
shock front followed by a variable stress field with an essentially elastic 


‘eehevier. The velocity of propagation of the shock front is 


i on the initial density ofthe material and E is the secant modulus of 


: ie are the normal vertical stress and strain immediate y 


) 


which and 


= 


ly 


Int terms of the modulus E of the he line BA, F Fig. x , €g Cal can be expressed | by 


“a “Waves and Shocks in Locking and Dissipative Media iby M. G, Salvadori, 
alak, Weidlinger, ‘Vol. No. EM 2 2, April, 1960, 


— 
Sas 
gg 
$ 
= — 


WAVES: 


_timet t its initial ‘position, 


it follows that 


By) virtue of Eq. 8, 12, 18, and 14 M4, ‘the field behind the shock wave can be ex 
pressed in the form 


= f(z ct) e(z + +e ct). 


— 
5 
. The equation of motion of the material behind the shock wave is given |! RR 
Compressive strain isdefined by 
a - Eliminating wand o in turn from Eq. 8 and 12 gives the one-dimensional wave © 

— 
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In a space t time plot, the locus of me shock raged is a curve as shown in Fig. 4. _ 


Consider er a general 
labeled 1 and 2 on the shock curve e located by charectertatics passing endl : 


the field —-, By adding Eq. 15° 5 and 16 for the the point 2, it is found i: 
wis 


since @- -c t) is 
field point | @, t) 


ly, considering point it is found that 7 


are obtained. 
the locus of the shock wave, Zg(t), is known, then. Og and are given 
_ by Eq. 6 and 7, and the remainder of the field can be found by Eq. 20 and a 
The solution is by establishing, Zs(t), fora given surface pressure 
Consider a field octet: which is onthe surface z = 0 such as point B, Fig. i: 
|For this point the stress o given by Eq. 20 is the applied pressure p at time t. 3 
By substituting the values for 01, 02 and dy, ug 6 and 7, an 


is the given ontieoe pressure history. The heights shown in : Fig. . 4 labeled 21 
and Zg and the times ty and 2 may be also iene to be functions of = 
The time ty is related to t by 


7 Differentiating ‘Eq. 23 with respect to t 


= 


obtained, in which the derivative is to to be interpreted to be interpreted as in 


f 
s | constant along the characteristic, it follows that for the a. 
_— and by substituting Eqs. 18 and 19 into Eqs. 15 and 16 
= 5 (ty +g) - (01 - (21) 
"eae 
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Transforming Eq. 24 into the form > 


and defining z! by 
and using Eq. 


is obtained. Similarly fro 


Substituting Eqs. 27 and 29 


in which pis the function p(t) andthe primes denote differentiation with respect — 
_ to the variable t. Eq. 32 can be integrated with respect to t to yield the final 
equation 

(33) 


= 


33 m 


EM3_ 
dt c+, a 
4 
(26) 
a 
4 
Ise of p(t). The boundar 
Eq. ty is known. — 


‘Sary to establish Zz; atsome early timet as a starting point. For this purpose, 


short time approximation is developed in | Appendix, I, 


‘NUMERICAL /RESULTS 


ov 


The magnitude of the front stress Os and particle velocity 

F peak values experienced at a given depth), can be presented for all bi-linear 
_materials in the domain 0,/€ cSE < «by a family of curves with the single 
- a parameter €e/€c¢ as shown in Fig. 5 and Fig. 6. In the general case in which 
Oc > 0, the peak stress and velocity propagate at the constant values 0, and 
Uge indefinitely after the stress at the fronthas decayed to this intensity. This 
behavior is due to the | elastic action of the medium below the critical stress. 
’ level and is indicated by a horizontal dotted line in Fig. 5 and 6, aaa aataal 


ing to some typical val value an and Us Of the medium, 


_ Numerical aaa es are plotted in terms of a non-dimensional front stress 
s/Po for a ‘decaying pressure load of the form 


in which to is the relaxation time of the pressure. The location of the wave 

front, where the peak stress is experienced, is measured ~—s the n non- dimen-— 


sional depth z/ ‘(coto) . in which 


is the initial shock velocity of a locking medium corresponding to the limit 


> ow, The _non- dimensional particle velocity is defined by u s/o 

in which cog €¢ is the initial particle velocity of a locking medium. 

_ In the limiting case €e/Ec + «, which is reached by letting €¢ >0, the — 

medium responds | elastically, which is shown by the constant peak stress and 
particle velocity at any depth. The other limit, corresponding to locking be- 7 

havior is reached when the parameter Ee 0, which corresponds the 


ated with any, value of the parameter €e c approach the curves s representing 
the locking medium, The depth, beyond which the peak stress and peack veloci- : 
_ty of bi-linear and 1 locking medium becomes nearly coincidental corresponds | to 7 

the: time a at which the first and principal rarefaction wave from the surface 
—* the shock front. ad the surface pulse decays in a time of the — 

tor the > depth o of coincidence Ze is is given 


This approximate coincidence of peak: values is of practical significance, since 


it demonstrates that stress attenuation beyond a certain, determinable depth, | 


i 
— 
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(34) 


LOCKING MEDIA 


9025 10” 


_—ATTENUA TION OF PEAK PAR TICLE VELOCITY 
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Sane, 


depends ‘principally o on the he single parameter te. Furthermore, in the case 
where the stress-strain diagram is a continuous curve, | the decay of the peak 
_ stress and particle velocity can be estimated by a bi- linear approximation of : 
the curve, which in turn reduces the problem to that of a single parameter, 
provided that the val values of of lie beyond th the depth 


_ The behavior of granular or porous materials such as sand, under confined 
compression is characterized by a non-linear strain hardening stress-strain 
curve, with large permanent strains on unloading. © if the unloading occurs on 
a steep linear path, wave propagation problems in such media require extensive 
z numerical work. Such materials, however, can be approximated by a bi-linear 
ve stress-strain diagram with resulting simplification in the mathematical and 
numerical treatment. This ‘simplification | provides means to estimate the in- 
fluence of various mechanical properties of the medium on the attenuation of 


i surface pressures, and of Particle velocities. _ These investigations show that 


various bi-linear materials becomes a indistinguishable, resulting 
further simplification of the problem. 


Development ‘Command, ‘USAF. 


ituting Eqs. 37, 38, and 39 Eq. 32 and equating of each 
sides" of the 1e equation 


10 
] 
al Lifes allenuallion iS Intivenced ine Character Of the Tress- Tall 
stigations sponsored by the Air Research 
4 APPENDIXI. 
surface pressure p(t), such as shown in Fig. 2 can be represented by 
values z; and may also be written in series orm 
24(t) = ky t + Ko t +. ‘ eee (3 
‘Subst 
i= 
0) 
— 


STRESS WAVES» 


function z1 (t) must also be such that at the appropriate time it is equalto 


_ the value of zg at an earlier time. may be written as 


3 which all symbols in brackets represent independent variables ofa function. = 


Substituting Eq. 38 and 39 into Eq. 42 
ky t2 + + 


Equating coefficients of eachp power of ti in n Ea. 3 37 after expanding» each term 


gives 


_ Eq. 40 and 44 may be solved for 


4 If the first two: yr in ‘Eq. 37 are considered adequate 1 up to a time t, the 
_ corresponding value of z1 is given by the first twoterms of Eq. 38. The value 
of ty is given by Eq. 23. These values form a starting point for the computa- 
tions | using Eq. 33. After z9 | is established, by Eq. 33 for a given t, , the value 
tg is found by Eq. 28. ‘The computations proceed with these values | of of 29 il 
te of t the previous § step as for the next 


ie: 
| 
Eqs. 4land45 maybe solved for ke and m9 considering ky and my known 


¢ = acoustic velocity; 

= secant modulus; 


of surface pressure; 


= constant; 


pressure; a ue 


= location of shock front: 


= density; and © 
stress, 


a In conducting numerical computations, IT IS CONVENnIent tO 
&g as Eq. 33 so that the numerical results for Z9(t) developed from Eq. 33 need [im 
| APPENDIX II. - NOTATION 
&g — 
Po: initial value of surface pressure; so » 
b 
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VISCOELASTIC WINKLER WITH INTERACTIONS 


The concept - foundation, consisting of. ‘the Winkler- 
foundation with | ‘shear interactions, is extended herein to the case of visco- 
elastic behavior. ‘The D. E. for. the vertical surface displacements w due to 
creep is formulated and solutions for several loading cases are worked out. 
In addition to the elastic constant two material parameters enter the analysis: 
The viscosity parameter p related to the shear deformations and na . visco- 
compressibility parameter of the vertical foundation elements. A procedure 
determine these parameters is suggested. The of this 

to the behavior of snow foundations is 


a such as compacted snow masses or certain frozen grounds show large creep 
deformations when subjected to loads of long duration. One of the main 
characteristics is that these deformations vary with time and nd by fa far exceed 

Several papers have appeared on the subject of beams or plates” ona 


. It was observed in the field as well as from experiments that materials | 


- Nisco-elastic foundation.2,3 In these papers the Winkler assumption was ex- 


- Note.—Discussion open until November 1, 1961. To extend the closing date one © 
month, a written request must be filed with the Executive Secretary, ASCE. This paper — 
is part of the copyrighted Journal of the Engineering Mechanics Division, Proceedings 
3 of the American Society of Civil Engineers, Vol. 87, No. EM 3, June, 1961, © 
1 Asst. Prof. of Aeronautics, New York Univ., New York, and Consultant to Snow Ice 
Permafrost Research Establishment, Corps of Engrs., U. S. Army. 
_ 2 *The Infinite Elastic Beam | ona Linear Viscoelastic Foundation,” by A. Mz Freu- =a 
denthal and H. G. Lorsch, Proc. Paper No. 1158, ASCE, Vol. 83 , No. EM 1, 
ee 3 “Flexible Surfaces on Viscoelastic Subgrades, . B. Cc. Hoskin and E. H. Lee, | 
Proceedings, ASCE, Vol. 85, No. EM 4, October, 1959. ee 


— 
: 
4 
INTRODUCTION AND STATEMENT 
| 
a 
| 
: 


June, 1961 


qm in 


ASE. TIME- SETTLEMENT GRAPHS 


elements nor like an isotropic continuum but somewhere in between. It was 

found namely that deformations: due to. outside loads are confined to the 
. vicinity of the loaded region. To bridge the gap between the two extreme 
1 foundation models, M. Hetényi,4 F. ASCE, suggested imbedding a continuous 


“Beams on Elastic M. Hetényi, The University of Press, 


tended by adding linear viscous elements to the elastic springs inparallelor | 
geries arrangements (Maxwell model, Kelvin model, or their combination). 
interaction between the partic ndation elements was assumed. 
. Experiments indicate that — 
behave like a Winkler-ty 
| = 
‘FIG, 1.—FOUNDATION MODEL ACCORDING TO HETENYI 
46 arts @ @w:! g ele ao 
2, FOOTING LOAD TEST ON A SNOW 
— 


assumed that the beam deforms in bending only. P. L. treated 
problem by assuming the existence of shear interactions between the spring 
elements. This approach is equivatent to the assumption of a beam or plate 
that rests on. a spring foundation and deforms only due to transverse shear. 7 
With this assumption the foundation is characterized by two constants; — 
representing the elastic foundation, the other the intensity of elastic shear 
~ interaction. From the point of view of the mechanism of deformation of com- 


34 FIG. 4. 4.— FORCES ACTING ON A SHEAR LAYER ELEMENT | 


-pacted snow masses, it seems that the - foundation is avery s suitable 


a Herein the concept of the Pasternak-foundation is extended to the case of 
, and the behavior of the foundation surface subjected 


7 to different load situations is studied as a function of time. Results of footing | 
: load tests performed on asnow base by the Snow Ice and Permafrost ma 0008) 


Establishment, Corps of Engineers, U. S. Army (SIPRE) in Greenland (1959) 
are shown in Fig. 2. It canbe seen that, after a relatively short initial period, 


5 “Osnovy Novogo Metodo Rascheta Rundamentov na Uprugom Osnovani pri Pomosh- 
chi Dvuch Koefficientov Posteli,” P. L. zd. Lit. po Stroitelstvu 


_ Architekture, Moscow, 


FOUNDATION, 
a q 
> 
a 


time displacement relationship exists. one prerequisite 
is that the foundation dimensions should be large in comparison to the dimen- 
sions of the loaded area.) This suggests the possibility of assuming Newton’ s 


ship. Because the elastic ‘deformations, that take e place immediately after | the 
load application, are very small as compared to the creep deformations, it - 
justified to assume a foundation as shown in Fig. 3. total deflection 


Ytotal ~ Welastic (1) 


Whereas Welastic is in the usual manner, the sections 
will deal with the deflections w due to creep. 
N eaten. .— The letter symbols adopted for use in this paper are defined 
whe: ar, in the illustrations or in the text, and are arranged 
FORMULATION OF THE PROBLEM 


consisting of incompressible vertical elements, “that ‘deforms by 
transverse shear only, and that rests on closely spaced linear viscous ele- 


m ‘and iso ropic then 


“surface | are slow, the inertia terms are 
The equilibrium equation 
° 


_ 


aw 
in which is the « downward velocity of the foundation surface. 


oa should be noted that two- constants of the foundation material | enter the — 
analysis: viscosity coefficient ‘related to. the shear r deformation of verti- 
cal foundation elements and nN, a visco- compressibility coefficient. Substi- 


-tuting 4, and 5 into Eq. 2 


— 
oA 
a 
| — 
i} cut out by the vertical surfaces x, x+dx, y and y+dy as shown in Fig, 4 _ iat 
— 
to Newton’s la 
4 i 
— 
— 
« 


if the load is time independent, p 
with respect to t it ete to 


\ox? 


Eq. 7 is the D. E. for the ential displacement of the foundation surface. 
- According to the assumptions made, the displacements in the x-y plane are 
zero. It should also be noted that no restrictions w were imposed ‘upon the the mag- ” 


nitude of w to obtain as a linear equation. 


‘setting the the stresses 
| 
as was done by E. Reissner” who iecindiaie the — of an elastic plate 
- resting on a visco-elastic foundation of the Kelvin type with shear interactions. 
In the following, solutions of Eq. 7 are evolved for a number of practical 
loading cases in order to facilitate proof of the applicability of the preceding © 
om for different foundation materials — a series of tests in the field, b 


x, t) and Eq. re- 
x2 


——— ‘the solution in the f ‘the form w= X(x) TW), Eq. 9 will b be satisfied if 


(11) 


wet deflection w is hey sian to vanish at Xx: = = 0, , SO that A= = 0, ane because 


7 * _ 6 “A Note on Deflection of Plates on a Viscoelastic Pe,” by E. Reissner, 
Journal of lied Mechanics, March, 1958. 


WINKLER FOUNDATION 
— = p(x, y), after differentiation of Eq. 6 _ 
= 
— SS 08 
wx, t) = (Cy + Cot) (Ae q 
q 
if 
q 


FIG. . SECTION THE RIGID STRIP | BASE 


g 
— 


) 


and 
Thus the surface of the medium will deform a according to an exponential func- 
tion (see Fig. 5) and the time deflection relation is linear and will be of the 7 
nature ‘shown in Fig. 6. The discontinuity in slope along the line of applica- 


tion of the’ “concentrated” line load is a consequence of the assumptions made | 
and will disappear once the sine — is seni wale one of finite width as done 


UNIFORMLY p IN THE INTERVAL - +b 

ee The solution for loading case of the u uniformly distributed load p in the 
interval -b = x = * is obtained from the previous result setting P= pdé and 


integrating from - -bto +b. The deflections” the Ic loaded region x 2>bare 


the deflections inside the region x<b 


FOUNDATION LOADED CENTRALLY BY A RIGID S sid 


section the ‘rigid strip base of a foundation loaded by the 
rigid strip is shown in Fig. 7. Under the strip, -b < x < +b, the deflection | WO 
is constant. ‘The | ‘corresponding vertical is, to Eq. 5, 


in 
The deflection suriace for x b is determined the continuity condi- 
that at time t, (x, , t) will pass through wo(b, t) for any A. Thus, 


% to 2, at x = bo 


follows that C = C ir in n Eq. 12 
= 


4 
EM3 WINKLER FOUNDATION 19 
4 ‘The constant C is determined from the condition of vertical equilibrium for ik 
(With: 
2 
| 
“4 
tm 


t length of strip) 


This is the velocity of a downward motion a the rigid strip of width 2b. 
if it is assumed that at = -0 the starts. 
, _ Because the deflection wo under the strip -b < x < +b is constant, no )shear — 
forces oc occur in this under the e strip is 


because e of ‘equilibrium (Eq. on two line loads 


per unit length of strip act along th the edges x = +b, as shown in Fig. 8. 
RIGID STRIP SUBJECTED TO A CENTRAL LOAD AND MOMENT | - 
_ Because of the linearity of the problem involved inthe case of a rigid strip 
subjected to a central load and moment the vertical displacements can be a 
‘ into a downward translatory motion and a rotation as shown in Fig. 9. Each 
_ part can be treated separately and then the displacements added. Peek 
A condition is that the total load ‘should not create tension in _ > contact 


It is assumed that the load P and the moment I M are time independent. ——. 
The case of the rigid strip exerting a central load was treated previously. 
The problem of the rigid strip subjected to a moment M will now be solved. 


_ With the notation as in 


Fa 
* Substituting in the vertical equilibrium equation (per uni 7 a 
oor 
&§ 
— 
= 

— 


WINKLER F FOUNDATION 


=. 
-- 


“ 


_-SUPERPOSITION OF EFFECTS DUE TO P P M M 


off 


g 
— 8.—PRESSURE DISTRIBUTION UNDER THE STRIP 
4 
— 10,—A SECTION THROUGH A RIGID STRIP AND BASE SUBJECTED 


4 for x >b in which angular velocity. = 
Moment  equilihriam ab about x=0 requires that 


The angle o of rotation ata time t is the 


M 


assumed that at t= 0 creep starts. 


_ FOUNDATION SU SUBJECTED TOA CONCENTRATED LOAD, P 7 
A fc foundation subjected to é a concentrated load. P is is an axially symmetrical 


The solution of Eq. is now 


Because w is expected to aoe zero at ay ‘it t follows in 0. 


deformations due to creep are zero, thus 0. "Letting 
w=Ct Ky (A r) 


The constant C is determined: the equilibrium condition. With 
nC Ky ( 


a 
Substituting Eqs. 2§ elds” 
rt=0 
: C yields 
4 
4 
= 
3 
&§ 
J J ardrdt=297 cs r)dr=27C 2 + + (39) 
7 
— 


The graphical representation of the deflection function for ‘different 
is shown in Fig. 11. Eq. 41 contains a atr= = 0. 


or 
‘FIG, 11,—DEFLECTION FUNCTION DUE TO A CONCENTRATED LOAD 
oe 
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deflection v will be finite, however, once the load is assumed over 
a finite area, as shown in the following section, 2 ees 
‘FOUNDATION SUBJECTED TO A UNIFORMLY LOAD 


INSIDE A CIRCULAR AREA | 


—_ For the case ofa foundation subjected to : a uniformly distributed load inside — 


a a circular area, consider an area A of the surface subjected to an arbitrary _ 
load distribution as shown in Fig. 12. 


The deflections at an an arbitrary point M- caused by a load 9 at 


=q dA =q (a, 


d 
) dp dg.... (43a) 
and que to ‘the entire distribu ted load 


al and the integration is to over the loaded region A. Assuming q 
° 4 isa constant distributed — a circular area of radius a with center at r = 0 
yields 


(r, t) = a‘ Ko (ar) fp p) dp 


and 1 applying fc formula 835. 3 3 given by | H. _B. Dwight” 


t) =q (A a) ) Ty a) Ky (? for r>a 
‘The deflection expression for. inside the: loaded ‘region 
‘Fig. 
t) 


~ Applying Eq. TT (see Appendix II) to the first integral expression and Eq. 


“Tables of Integrals and Other Mathematical ‘Data, ” by 8. , Revised a 


iim 
Bet 
Using Eq. 77 (see A 
2 
| 
i 
is (see ig 
4 
pf. .(48) 
718 
(49) 


Using formulas 835.3 and 835.5 given ste ontitw and with the identity 


1- a) K, ( 


‘FOUNDATION LOADED CENTRALLY BY A CIRCULAR RIGI RIGID PLATE 
section through the rigid plate of a foundation loaded centrally 


by. a circular rigid plate at atime t i is as replacing b D by rad-— 


ius a and x byr. By : an argument similar to that previously given under the 


section headed “Foundation Loaded Centrally by a Rigid Strip,” denoting for 
any t the corresponding vertical displacement of the iad by wot), the de- 


flections in the region | r Saare 


We = Wo 
> 


| 
4 


“Vertical equilibrium requires that 
27 a 27 
J ar rdr de = =0 
Substituting the im for q into Eq. 58. and integrating yields 
+ 
+2 Ky Aa 
is of rigid plate. ‘The corresponding vertical dis- 


if it is assumed that at t = 0 ee 
The corresponding vertical in the contact area a of 
plate and base is similar in nature tothat described inthe section dealing — 
i foundation loaded centrally by a rigid plate and shown in Fig. 8 with the | uni- 


ay 


dt / 


CIRCULAR RIGID PLATE TOA A 
-- CENTRAL L LOAD AND MOMENT 


: The ‘procedure for a circular rigid plate subjected to a central load and 


and Moment” and illustrated in Fig. Here ethec condition n that no ten- 
— 


s 


which a; a is the radius o "Here also we assume that the load 
P and moment M are time independent. 


Because the plate loaded has been treated previously it is 


"moment as shown in Fig. 14. 
‘The in the region rs safora known are 


; 
A 
4 q 
4 
along the circumference of the 
| 
= = F 
4 


and the distribution is 


which: Wo = da, is the e angular velocity of the rotation. Now the 


tion surface outside the rigid plate ( (r, > 2a) ist is not axially ‘symmetric. For its: For its 


= 


rmination Eq. must be used in the form 


shape of the deflection surface Wj at a suggests that the detections 
a 2 r will be of the form 


= 


2 
q 
_ 
| 
2 


‘The of is then 


We is to be zero at infinity and because at t de- 


formations due to creep are zero setting C2 B=C 
We = Ct Ky r) 


‘The ec constant Cis determined : from the condition that for any tat T=a, We = Wie 


Thus 
os 


ne s (Fig. 14) | requires 


‘Equilibrium of mom 


0) rar do (r cos 


fe) 
‘Eqs. 63 and into into Eq. 70 and yields 


a)K, (xa) +4 Se (A a) 


4 
Ky 


Ky (0a) +4 Kp(ra)] 


r dr =0. .(70) 


The angle o of. rotation Wo at a time t is then 


aK, (Aa) Mt 


at a +4 Kp (0 


‘The determination of | m and n can proceed along lines suggested by Paster- 
nak for the elastic case. For the case treated here, namely time ‘dependent 


deformations, plate of radius a is tobe subjected to a constant 


load eccentrically placed, with ai an eccentricity e<*. From the time deflection 
t versus t are ‘The: 


4 Wo v 
dw, 


rar 

+ BK,(\r)| cos @..... (66) 
i 
ae 
if it is assumed that at t : 
‘ 
— 
rif 
| 
at 


veloped theory for the foundation material under consideration. = | 


al In this connection it should be noted that the plate diameter should not have 7 
~ any influence upon the determined material constants. Experiments should 
also establish the maximum contact pressure for which the time-deflection 

relation is still (nearly) linear, because this is one of the conditions on which — 


‘This be first to. prove the applicability of the de- 


ACKNOWLEDGMENTS 


the’ writer wishes to thank W. K. Boyd, F. A 
_ Branch, SIPRE for helpful comments and A. Assur, " Applied Research Branch, 


SIPRE for a critical review of the manuscri 


I.—NOTATION 


The following symbols adopted for use in this paper a: are li listed h here f for 
‘convenience of reference and for the aid of discussers: 


A, B, Cs, = integration constants; 


radius of loaded d area or circular rigid plate; 


ae = width of loaded strip or loaded rigid sees 
= excentricity of load; 


external line load; 


= distributed load; 
= vertical pressure of viscous: elements shear layer; 


= vertical pressure outside the loaded region; a 


= vertical pressure inside the loaded region; 


constant vertical pressure in the region of rigid plate; 


= line reaction along edge of 
= polar coordinates in the surface of foundation; 


s deflection foundation s surface perpendicular to x,y y pane: 


b; 


We = deflection the loaded 


a EM 3 
- 
4 
| 
4 
= 


= deflection inside the loaded region; 


= constai t deflectio under rigid late; 


cartesian n coordinates ir in the surface of foundation; 


= viscous foundation an d wre 
Migs ¢ 


E K,(AR) - FUNCTION 


and because of symmet ry 


R)= pI, 


™ 


Ky (Ar) p) 


8 “Theory of Bessel Functions,” by G.N. Watson, Cambridge Univ. Press, New York, 
1945, p. 361, Eq. 8. 


4 
a q : iq 
— 
= 2 and 
‘ 
(AR) doé=27 K,Qp)I,Qr)  forr< (78) 
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ULTRASONIC EVALUATION OF \F REINFORCED PLASTICS 


‘By Albert G. H. Dietz,1 M. 


a Reinforced plastics are becoming engineering materials of increasing im- 
_ portance. A critical problem is todetermine by nondestructive means both the | 


quality and composition of the material in finished structures. 
a Bursts of ultrasound ranging from 1.5 megacycles to 2.1 megacycles per sec 
and of afew microseconds duration were transmitted through plain unreinforced 
"polyester and epoxy castings and through specimens of laminates with 25 plies 
of 185 glass fiber fabric per 1/2-in. thickness, several thicknesses of — 
ranging from 1/4 in. to | 1 in. thick, laminates made with varying amounts of en- 
trained air, and laminates" with different finishing treatments applied to the 
» Velocity. of transmission of the ultrasonic pulses proved | to be directly pro- 
portional to the ratio > of glass to resin; the higher | the glass content the higher 7 
the velocity. Attenuation of the signal passing through the specimen was pro- 
‘portional t« to the amount of entrained air, that is, , to the soundness of the spec- 
imen; the higher the air content the greater the loss in energy of the signal. 
Finish on the fabric had no consistent effect. 


Some tests were made with close coupling by silicone oil and the transducers 


- pressed tightly against the specimen, others were made by transmitting the . 
pulses throughan thickness (1/2 in.) of water. Both methods gave 


7 _ Note.—Discussion open until November, 1961. To extend the closing date one month, 
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vely concluded that the method shows promise of 
destructive test for quality a and composition of reinforced plastics. 


INTRODUCTION 


"amounts £ for structural applications. In use: materials the glass fibers im- 
part most of the strength and the resin acts as both binder and support for the 
fibers . The structural properties therefore depend on the ratio of glass fiber ; 
and resin. The soundness of a reinforced plastic part depends to a large extent 
on the care taken wn during fabrication to avoid porosity, because porosity ~~ 
permit» penetration by water or other fluids and this, in turn, may lead to de- 

terioration, 
_ Visual examination of a reinforced plastic part generally | cannot reveal i 
‘ratio of glass fiber and resin present, nor is such examination a reliable indi- 
-eator of the quality of the part. Some nondestructive means of ee 
composition and quality i: is consequently highly desirable. 
a _ Of the various nondestructive means of investigating materials, it ne 
- that propagation of ultrasound is one of the most promising for the nie 
of reinforced plastics. Because the velocity of transmission of a sonic all 
through glass: is higher than through the polymeric. binder, increasing glass 
content ina specimen should result in increasing velocity of transmission. A 
“pulse passing througha porous specimen containing air bubbles should be scat- 
tered internally andits energy considerably attenuatedas compared with a pulse’ 
passing through a nonporous specimen. Measurements of time of transmission © 

— andof attenuation of pulses of ultrasound s should therefore be bei indicators of com-— 7 

_ To test this hypothesis, a series of e: experiments was conducted under the 

_ sponsorship of of the Naval Ordnance | Laboratory, White ee . All spec- 

imens were prepared by at agency. 


in Fig. 1, area | rectangular pulse generator that feeds pulses of approximately _ 
5 microsec duration into a barium titanate transducer on one side of the spec- 

_imen, a similar pickup transducer on the other side, a pulse receiver with | 

attenuator, and an oscilloscope. ‘The shift of the received pulse relative to the 

_ transmitted pulse on the time axis of the oscilloscope is a measure of the time 
of transmission, and the change » in the amplitude of the pulse is a measure | of . & 
attenuation of the signal as it passes through the specimen. ieaigiiiaieeiitl 7 
Two types of coupling were used. Inthe first series of tests the transducers 

; were pressed directly against the sides of the specimens in a film of silicone 
oil to provide close coupling. In the second series the specimens were im- 
-mersed in a water cell and the gesagt also immersed in the water cell, 7 

_ were held approximately 1/2 in. away from the ‘sides of the specimens. The” 


_ water-immersion technique was tried to see if good results could be obtained 


“ported®,3 and will not be described in detail here. . The essential features, shown 


a 2 “A New Tool for Determination of Stage of Polymerization of Thermosetting Poly- | 
mers,” by G. A. emia and E. A. Hauser, Journal of Polymer Science, Vol. VIII, No. 6, 

3 “Cure of Phenol- For 1yde 


> 

7 

> 

4 
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EM "REINFORCED PLASTICS. 33 
with transducers some distance away trom the test object. In any practical 
: test this would be much simpler than attempting close coupling with silicone oil. = 
_ Inthe silicone series panel specimens (described subsequently) were first 
_ tested dry as received. Panels were then cut in halves, and 1/2 of each ganel 
immersedin water, heated to for 1 hr, to cool gradually, and 
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maintainedat 24°C, fora total of 48hr. The objective of this | Series wai was to see 
if moisture pickup would have a detectable effect onthe transmission of ultra- c 
sound. The other halves of the panels were reserved for tests not pertinent to 
those reported h herein. 
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= 
| 
| 
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ae TABLE 1,—ATT ATTENUATION AND VELOCITY OF ULTRASONIC 
CLOSE COUPLING AND WATER IMMERSION 


Variable 


Mean deviation, 


Content 
Glass Content 
Glass Content 
18 Plies 18 Thickness (1/4 in.) 
68 Plies 18 Thickness | 
35 Plies 182 136° 
35 Plies Maximum ‘Air 
35 Plies 182 | Volan | Finish = 


Polyester Casting 


21x 10® cycles per sec 


Standard iation x 100 = of variation, 


Standard Deviation 


Two sets of fr equencies" were employed. In the first, a fixed frequency of 
2. 1x108 cycles per sec was usedon all specimens. In the second, the frequency 
was varied from 1. .5x106 to 2. 1x10 cycles per sec with certain of the spec- 
imens to see if any ‘frequency effects might be noted. This. might be particu-— 
larly noticeable if under some combinations of frequencies and compositions 
internal standing waves were generated. 
Test Specimens. —Specimens were prepa aredto make possible a 
of the effect of the following factors on ultrasonic transmission: s => 


Ratio of resin and glass fabric, 


2. Thickness, 


Type of finish employed on the glass" fabric. . (Because polyester resins 


5 2 not inherently bond well to glass fibers as manufactured, various ty pes sof - 


finishes are applied tot the glass to improve the bond. Po 
Glass fabric was type 182. 


_ Finishes employed were those known in the trade as 136, Volan, and 112. a 


- Except for the series of panels to test the effect of thickness, all panels were 


_ For purposes of comparison, 2 panels of cast resin without reinforcement > 
_ were included. One consisted of the polyester mix employed in the test —- 


imenand contained 80% ofa hard rigid formulation, 15% ofa softer more flex- 
ible type, and 5% ‘styrene. (Rohm and Hass Paraplex P43 80%, Paraplex P13 — 


= and Styrene 5%.) ‘The second cast resin panel was a standard epoxy. 


(Shell Epon 828) 


tm 
‘Panel No, | Glass Finish 
: 
2 190 | 
140 | 4.47 
8 15.5 | 4,44 
&§ | 1.61 
12 3.9 4,99 
130 
: 
f 
| 


(dbperemy "Velocity (105 en (105 em per sec) ater 


Wet 


deviation, 
iat cm per sec) 


nwo 
oww 


nwo 


NNNN NN WW 


Panel: characteristics are summarized in Table ie 

Panels 1-5 had glass fabric varying from 25 plies to 45 plies, 3, consequently 
the ratio of fabric to resin varied but otherwise these panels were identical. 
a Great care was taken to minimize air entrapment. . It proved impractical to 

make good panels with fewer or more plies. 
Panels 6 and 7 were 1/4 in. and 1 in. thick. V With Panel 3, they constituted 
_ In panel 9 as much air as possible was deliberately er entrapped; in Panel 8a a 
7 - medium amount of air was added. Panels 8, 9, and 3 (minimum air) constitute 
a series onair entrapment. It was not possible to _— the actual amount 


_ Panel 3 with 136 finish, Panel ss with Volan, and Panel 11 with 112 finish — 
constitute a series on finishes 
¥ Panels 12 and 13 were cast resins. Panel 12 was the same resin as used in 7 
Panels 1 through 11. , Panels ‘12 and 13 afforded comparison between two resins. 


Results. 
1. 


‘ag ASTICS 
a 
Mean 
98 (2,16 1.3 | 31 
16.09 1,99 16.4 2.83 
og 18,52 5.51 2.95 16.9 2.63 
15,64, 2.85 | 2,97 = 
4 
g 
; 
4 Dry and Wet Panels.—ResultsaresummarizedinTable 

is - and 10. Columns 6 and 8 require a little explanation. | te 
Standard deviations were computed in the usual manner. 
a 
ie. ; in whichS isthe standard deviation; @ represents the individual value; @ denotes ae 
a the mean value; and n is the number of observations. To make comparisons a 


June, 1961 
easier, the standard deviations are given as ‘poecentage of the mean. | This sis 


“known as the coefficient of variation.” 


The salient points are the following: a 
- Panels 1 to 5 showthat with increasing fabric content there is a decrease 


in inattenuation andan increase in velocity. As shown in Fig. 2, the velocity pro- 
"gression from 45 plies down to about 30 plies is quite linear, but then the dry a 
value deviates to a curve that canbe extrapolated tothe dry value for zero plies 
% (specimen 12, Table 1, cast polyester) with its value of 2.54 x 10° cm per sec. — 
The wet velocity value for 25 plies deviates from the otherwise regular pro- - 
eect . Because no wet value was obtained for the cast polyester block, the -" 
= extrapolated curve for the wet condition is only estimated, and shown dotted i in 
2. Fig. 3 shows a fairly regular progression of attenuation values, maxi- 
mum for 25 pliesand minimum for 45 plies, except for 30 plies where the val- ‘J 
ues deviate considerably. It is to be noted, moreover, that the attenuation for : 


would seem to indicate thata maximum is reached somewhere above zero ane 


straight lines than do the values, the attenuation values. of the wet 
are consistently higher than those of the dry, and the velocity values are con-— 
‘sistently lower, 
as 2. Panels 3, 6, and 7 when dry show a regular pr progression of attenuation — 
- with thickness. This is brought out in Fig. 4, in which the attenuations are — 
versus thickness and extrapolated to zero thickness. The line passes’ a 
oa through: zero attenuation, | that is, through the origin. . The wet values do not lie 
closely on a straight line, but the best average wet line lies above the dry line, 
that is, all attenuation values: are higher for the wet material than for the dry. 
oy: Velocity values are. fairly close, indicating that for the given ratio of resin 
and glass, held reasonably — constant in the three thicknesses, the velocity is - - 


reasonably constant. 

a. 3. Panels 3, , 8, and 9 show that with increasing air entrapment there is an 
increasing attenuation, as is to be expected from the high attenuation of air. 

_ Scatter in the dry results is also greater, as indicated in the greater values of a” 
standard deviation for greater air entrapment, but scatter shows no consistent — 
trendin wet values. Velocity values remain fairly constant, although lower for | 
wet thanfor dry materials, indicating that for the given ratio of resin to fabric : 
the velocity is not affected. 


4, Panels 3, 10, and show that different finishes affect the attenuation 


than dry. Volan developed the greatest attenuation, and 136 the least, although ~ 
it did not reed greatly, when oo 112, which in turn had the smallest 


o or dry, for agiven ratio of fabric to. resin although wet velocities are lower 


_ deviation than 136 or 112, but all three showed greater standard deviation ~se 
: wet than dry. No definite conclusions can be reached as to the effect of finishes. — 
5. Panels 12 and 13 show that the epoxy and the polyester mix develop ap- — 

_ proximately the same velocities when dry, but attenuation and standard devi- 

_ 4 “Statistical Methods in Research and and Production,” by on 3 Hafner Publishing 
Co., New York, Ed. » Ps 
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standard deviation than Volan. Volan and112 showed greater increases in 
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REINFORCED 


- ation of the eam mix are markedly lower than the epoxy. No wet tests 


were made, 
In the attenuation values are higher for the 

than for the dry. Velocity values are lower i in every instance but one ( Panel 7) 7 
in which there is no essential change. In Panel 4, likewise, although there is 
a drop, it is not significant. Velocity changes, ins general, are smaller than 


; Water Coupling, Dry Panel. .—Results are summarized inT Table 1, col- 

umns— and 12. These values were obtained | by ‘subtracting from the total 

elapsed time the time needed for the pulse to travel through the water (from 
= values of velocity of sound in water) and by similarly correcting for the | 
known values of attenuation in water. 

_ One unavoidable source of error was present . The previous tests were made > 
7 with Silicone oil coupling. . This oil could not be completely removed even by : 
a: intensive | cleaning. - Consequently, perfect wetting of the surface could not be 

- obtained. _ Imperfect wetting would be expected to lead to higher attenuation 


than would occur with good wetting. 


Resultsare summarized in Table 1 and in Figs. 5, 6, and 7. compari 


son, corresponding results obtained with oil couplingare included. Fig. 5 gives | 
the velocity values of 25- ply to 45-ply panels and Fig. . 6 gives the correspond- 
ing attenuation values. Fig. 7 7 gives attenuation . values for 1/4-in., » 1A- -in., and 


salient results are | the 


Som "Velocity values | in dailies coupling run slightly lower than in oil | coupling. 
‘This difference may be real or apparent. It has already been noted that a thin © 
residual silicone « oil film made perfect water coupling impossible. it may b be 
that this influenced all readings, including time of propagation 1 of the: pulse, and 
therefore led to apparently reduced velocity values.  — 
Although the water- coupling velocities are slightly lower than the oil- 
coupling velocities, ‘they are consistent and show the same trend as the oil- a 
_ coupling values. As shown in Fig. 5, the water-coupling velocity values are 
i more consistent than the oil, and the anomaly noted in the 25-ply oil- — 
velocity does not appear in the water- -coupling value. 
_ In the various triplets of corresponding panels, that is, Panels 3, 6, and se, | 
= 8, and 9; and 3, 10, and 11 velocity trends in oil and water are similar, and a | 
deviations are no greater than the experimental error, 
_ 3. Attenuation values in water coupling generally run higher than in oil. It 
“has already been pointed out that this effect may be more apparent than real : 
‘because of imperfect wetting of the silicone oil-contaminated surfaces of the 
“specimens. The effect m may be realif aging has caused degradation of the resin. 
‘These panels had aged approximately a year between the time the first tests 
were run and this series begun. Attenuation is particularly sensitive to the 
degree e of ct cure of resin as brought elsewhere.: 2 Degradation} about 


asimilar realeffect. 


4, Although water-coupling attenuation values are higher than 


‘they are consistent in themselves as brought out in Figs. 6 and 7 and in Table 
1. In Fig. 6 the 30-ply anomaly previously noted still exists but is less marked 
in water- coupling than in oil. The factthat the 30-ply attenuations arethe same © 
for oil and water is probably a coincidence o1 or an indication that the surfaces 
ofthis panel were somehow exceptionally clean of oil. Identical or nearly identi- — 
cal values are also found in some other panels, such as 7, 9, and 10. The ex-— 


a 
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- TABLE 2.—AT TENUATION AND VELOCITY OF ULTRASONIC 


n (db/en m) 


Second Series 


.7 


Glass ‘Content 
vin 


Glass Content 


| p-43 Mix | Casting | 


y= 1,5 x 106 to 2 1x10 cycles persec 


_ perimental errors, larger for attenuation than for velocity readings, can sleo. 

account for apparent coincidences or deviations from the general trend. 


me... 5. In Fig. 7 ‘the water- ~coupling line lies above the oil- ~coupling line. In 


. nat 3, , 10, and 11, the attenuation trends in water and oil are thiecd simi- 
lar, although apparent inconsistencies appear. Most of these are within ex 
perimental error. Thus, in 3, 8, and 9, oil-coupling attenuation for 9 is higher © 

than water- coupling, but the difference is smaller than experimental error of 
. The same holds true of the apparently identi- 


_ 3. Varying Frequency.—This series was run with water coupling to see if 
varying frequency would have a noticeable effect on the results. - Among t these 
_ was the possibility that standing waves might occur at certain frequencies and 
“that these would affect the results. Tests were run only on Panels 1 to 5 and 
_ Results are presented in Figs. 8, 9, and 10 _ in Table 2. For comparison, | 
the results of the dry tests inthe first series of tests are also includedin Table 
- 2in the first two columns headed “First Series, Dry.’ ” Fig. . 8 shows graphs | of 


foe frequency. Fig. 9 gives attenuation versus frequency. From 


the velocity values of Fig. 8 the numbers of wavelengths in the thicknesses of | 
the various panels were computed, and in Fig. 10 the attenuation values are ~ 


- replottedagainst the numbers of wavelengths in the thickness. . The higher the - 
, velocity the smallerthe number of wavelengths in the thickness for a given fre- 
quency; consequently, the high glass- ‘content panels, that show the gromtest 
locity, are shifted to the left in Fig. 10. 


_ The salient results are the following: _ 


_ eshte are e relatively ur unaffected by f es as shown in Fig. 8 and 
‘Table 2. increases S with increasing glass content at all frequencies. 


= 
1 | 25 Plies 182-136 | Glass Content 13.3 
ies = ‘ass Conter 7 7 
| 35 Plies 182-136 | Glass Content] 8.6 | 9.0 if} 123 136 
| 40 Plies 182-136 8. | 9.8 | 10.6 
loss at the interface for panels of this glass-resin ratio. This reinforces the vee 
- suspicion that the increase in water-coupling attenuation values over the oil en 
values is at least largely apparent rather than real. 
| 


“PULSES IN GLASS FABRIC POLYSTER PANELS* 


Velocity (105 cm/sec) 


Second Series” Ser ies, 


w 


oun 


aie 
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The values obtained in this series of tests are somewhat lower than those fe) 


tained in the first series, but are consistent inthemselves. _ a Hay 


kh Attenuation increases with frequency in all panels and is _ marked 
in the panels having minimum glass content. The cast resin panel shows the 
7 least attenuation and the least increase with increasing frequency. For the one 


frequency common to this series and the earlier series (2.1 megacycles per 
sec), the attenuation values are somewhat higher in the second s series than: in 

first, but the results are ‘consistent in themselves. 
_ The 25-ply panel shows more erratic results than the others, especially at 
is “the lowest frequency. Experimental errors, greatest when making attenuation 


_ in the first series may again have become noticeable here. The 30- ale panel, 
that deviated in the first series, shows almost the most consistent results sin 


3. If standing waves shits: Sebi to ‘marked re reinforcement of signal at 
= wave and extinction at quarter-wave and three-quarter intervals, their 
a should be brought out in Fig. 10 in which attenuation is plotted versus» 
the number of wavelengths in the thickness of the material. _ Although some of 
_ plots, notably for the 25-ply and 40-ply panels, show some up and down 
_ variations in attenuation, others do not, and the variations that do exist do not | 
appearto t be consistent with half-wave and quarter- wave or wave 


Be 


“sonic pulses passed through glass fiber-reinforced polyester panels offer a 
feasible nondestructive method of determining the proportions of glass fiber 


a 
‘First 
veries 
__Dry 
$i | 227135 1 129 146 | 22 j 28 220 2. 
| 2.63 | 2,63 2.63 2.63 | 2.64 | 2.64 
| 2.79 | 2.81 2.8 2, 
11.8 | (2.99 | 2.99 | 2.97 | 2.98 | 2.97 | 
6.3 2,39 | 2,39 | 2.37 | 2.39 | 
4 
not appear to be possible to correlate the variations with the fabric-resin struc- 
— rariation with wavelength seems to occur in the a 
4 
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F ‘1G. _ULTRASONIC ATTENUATION V VERSU -ERSUS 
-_ NUMBER OF WAVE LENGTHS, WATER | 


— 
= 


and te resin, and of deterniaiog | the quality of the laminate, v when where to 
corresponding standard panels of known composition ¢ and quality. msn 
2. Coupling may be satisfactorily accomplished | through water layers at 


‘tent 1/2 in. in n thickness, and does not require the use of thin oil films such | 


= 3. _ These tests unearthed no evidence that, in the range of frequencies em- ; 
ployed, the effect of standing internal waves is of sufficient magnitude to affect 
oe overall trend with respect to velocity and attenuation as affected by compo- 

_ Sition, quality, and varying frequencies. ” In this range, , velocity is independent | 
frequency, attenuation is ‘proportional to frequency. 
__ Betees the method can be said to be fully established more tests with differ- 
ent. compositions, , degrees of cure of resin, and with more e sets 3 of specimens 


“must be made. The procedure appears to be 
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WHITE NOISE REPRESENTATION OF EARTHQUAKES* 


Closure by G. N. 
— a e approximate derivation of the response spectrum 
_ _ presented by D. E. Hudson, G. W. Housner, and T. K. Caughey is interesting, — 
but the nature of their averaging process is not clear. The equation ee - 
shows an average oscillation increasing asymptotically with time to. In the 
_ ‘paper by Hudson, mentioned in their discussion, the time, to, is the time dura-— 
tion of the random pulses and therefore should be 25 sec instead of 6.5 oo for 
the case inquestion here, 
‘Their sentence t that states that the experimental points appear a bit — 


at the low partes end cannot be — ™ comparison with an approximate 


spectral the input isa 1 good one. increasing the energy 
towards the low period oun the the correspondence between the 


desired. However, as the discussers point out, further statistical ¢ comparisons 

: are desirable. Unfortunately, more information than is available is needed to © 
wed substantiate the situation. An indication of the frequency di distribution of 

— largest values 1 resulting from a series of bursts of white noise excitation is : 

_ shown in a further paper. 15 These results show a relatively small dispersion 7 

largest values. The chances of a value greater than twice the average value 
appear ‘negligible. In fact, approximately 90% of the largest values lie within. 


+ 20% of their mean. This is the order of the distribution of actual earthquake 
‘Spectra abo about their mean value 


“I April 1960, by G. N. Bycroft (Proc. Paper 2434). “rie PRE 
14 Engrg. Physicist, Dept. of Physics, Stanford Research Inst., Menlo Park, Calif. 
a 15 “Effect of Stiffness Taper in Aseismic Design,” by G. N. Bycroft, Bulletin, Seis. 


“Soc. of Amer., Vol. 50, No. 4, October, 1960. 
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BEHAVIOR OF BUCKLED RECTANGULAR ‘PLATES? 


Closure by Manuel Stein 


__ MANUEL STEIN. 21th discussion of the accuracy of the results p pre- a 


appreciably different results from the "subject It 
gratifying to know that an accurate analysis such as that of the discussers — 
gives essentially the same results. 


@ April 1960, = — Stein (Proc. 3445). 
21 Aero Research Engr., Natl. Aero and Space Admin., 
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BEHAVIOR OF VISCOELASTIC BENDING 


GEORGE 


‘The suggestion of using integral operators because of their close sainiioe 

tion with the creep function is a good one, Although the concept of solving an : 

elastic plate subjected to a suitable time-dependent loadl4 is a convenient 


one, it was not ado ted | inthis case. 


e* 


on a June 1960, by George E. Mase (Proc. Paper 3498). re 
_ 13 Assoc. Prof. of Applied Mechanics, Michigan State Univ. ., East Lansing, Mich. 7 

14 “Variational and Lagrangian Methods in Viscoelasticity,” by M.A. Biot, 1.U.T.A.M., 


a 
_Closure by George Mase 
oundary conditions would cer- 
&g - Tain y have added toward the completeness of the proportional loading solution. _ P, 
a : The expression in parenthesis to the right of Eq. 26 was intended to denote 7 ag 
> 
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‘MELVIN L. BARON, 11 M. ASCE and MARIO G. ‘SALVADORI,!2 F, ASCE.— 
The restraint method proposed by O. C. Zienkiewicz is particularly interesting 7 
whenever solution of the statically indeterminate problem of Stage 3 is already © 
_ available in the literature or whenever it appears to be easily computable. The — 
Boley procedure used by the writers would presumably be more advantageous _ 


whenever the solution of the statically indeterminate problem is cumbersome. 
‘The stress distribution in Fig. 5 refers to the crown of the arch, 
The writers are well aware of the fact that, as Zienkiewicz notes, the 


Duhamel integral presents a powerful tool allowing a solution of numerous 


Correction. _—The authors wish to correct a in | Eq. 30. should 


Ata 
June 1960, by Melvin L. and Mario G. (Proc. Paper 2516), 
11 Chf. Engr., Paul Weidlinger, Cons, Engrs., New Fork, N. Y. Y. and Prof. of 


_ 12 Assoc., Paul Weidlinger, Cons., Engrs. Nev New York, N. Y. and Prof. of Civ. Engrg. 


Columbia Univ., .» New York, ‘N. 


STRESSES DUE TO THERMAL GRADIENTS IN REACTOR SHIELDINGS® z= 
Closure by Melvin L, Baron and Mario G. Salvadori 2 
— 
a 
i 
‘= 


— 
‘ 
4 
' 


___ ERNEST T. SELIG,13 A, M. ASCE., ‘KEITH E McKEE, 14 M. ASCE, and 
_ EBEN VEY, 15 _F. ASCE,.—Regarding | the comments by Rongved, the 
fully recognize the complexity of the ‘soil-structure interaction problem 
buried structures. a nation’s national defense requires the design of such 
‘structures today and, hence, an attempt at solving the problem can not “be e 
deferred until well accepted theories of various specific buried structures" 
have been developed.” Furthermore such theories are most likely to evolve 
from a more general approach © that is gradually refined as details are con- 

_ The theory reported in the } paper was derived for the case of plane- sided 

a ‘structures, although extensions of this approach to arches, domes, and so 
forth, have been tried with some success . This theory represents an attempt 
to present a ‘rational mathematical model that explains arching of the soil 
over the structure. The soil was not assumed to act as a thick slab, but rather 
the equilibrium of the mass of soil in the zone of failure was evaluated to © 

; ‘determine the forces transmitted to the yielding structure under conditions of 


impending failure. In a sense, using this approach the soil has been treated © ; 


asa series of arches. To the extent that this has been done, the paper presents a 
a mathematical statement somewhat analogous to the shell theory outlined by ; 
-Rongved and having similar “great advantages” and disadvantages. 
A number of two-dimensional tests performed by the authors to aid in 
evaluating the theory” have substantiated this approach. Fig. 12 shows | the 
failure planes that formed when a rigid roof panel and a flexible wall panel | 
collapsed under a uniform surface surcharge. Fig. 13 compares the —- 
mental values of uniform, infinite duration overpressure Pg (0) required to 
collapse the roof panel at various depths, H with the theory. ———_” 
Rongved states that “the dome and shaped underground structures 
7 induce far more carrying capacity in the overlaying soil than comparable — 
flat- roofed structures.” It is difficult | to see how Rongved arrives at this. 
The concept in the “imperfect ditch” theory for tunnels under high fills — 
provides for sufficient deflection over the top of the tunnel to ensure the full 


ben yenefit of the load carrying capacity soil, It seems apparent | 


a August 1960, , by E. T. Selig, K. E. McKee and E. . Vey —_— ‘Paper 2575). | 


13 Assoc. Research Engr., Soil Mechanics Sec., Armour Research Foundation of 


14 Research Engr., Structural! Mechanics Sect., Armour Research Foundation of 
am: 15 Prof. of Civ. Engrg., Illinois Inst. of Tech. . and Supervisor of Soil Mechanics. Sect., 


Ar mour Research Foundation of Mlinois Inst. of Tech., Chicago, Wl. —™ all 
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Flexible Wall Panel 


b) Rigid Roof Panel 


FIG. 12.—OBSERVED FAILURE OF DRY SAND 
_ ‘SURFACE SURCHARGE 


_ 

UNIFORM 


DISCUSSION 

this ) provision a a flat ‘root w would undergo more vertical deflection than a dom ~~ 
“ hence, would permit more load to be carried by the soil. ‘play, 

Kovac and Frankian correctly | understand that the structural elements are - 
“assumed to have failed and that the soil strength is subsequently determined, | - 
The total resistance is then the ultimate strength of structural element 
plus the associated soil strength. There was no attempt made in the paper to — 
apply the to static problems eh they sh should 


DEPTH OF ‘BURIAL, H, (IN). 


FIG. 13.—COMPARISON OF THEORY WITH TEST R RESULTS 


> FOR RIGID PANEL ROOF 
It is apparent, however, that a buried structure does not have to support the 
full weight of soil above it and this fact is recognized in various ways in nor- 
, mal design practice. The ultimate strength of the structural elements were 7 
- therefore selected to represent what the writers felt was reasonable for — 
normal conditions. Hence, the ‘meaning of the symbol used for ultimate ae 
is the full resistance of the structural element. Big: ante 


sistance, Pm; 


“comes” ‘greater — “rate. appears to become infinite at 


| 
| 
7 
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finite “depth. “The in the paper exhibits this behavior as 
the _ experimental results as shown in Fig. 13. With regard to the soil pro-_ 
_perties given in Table 1, _ these were provided as stated in the paper for use 
if “other information is not available.” ” The comment that this was necessary 
“so that Eq. 21 would have significance” is incorrect, because for o= 0, as 


can easily be observed, Eq. 21 holds simply becomes 
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ARCH DAM ANALYSIS WITH AN. ANALOG 


by A. Hill 


RAYMOND A. HILL, 22 F. ASCE. —No attempt will be made by the writer 
to criticize or amplify this excellent paper on the application of eos 
analog computers to arch dam analysis, because the paper describes clearly — 
_the advantages and limitations of various methods of analysis. The writer _ 
wishes to point out only that preliminary designs must be carried out in 

more detail than would be the case if a digital used to 
analyze the stresses by the trial load method. 

_ On the other hand, it was found in the analysis of Blue Ridge Dam, — 
scribed by the author, ‘that variations in loading, changes in boundary condi 
tions, 
analog computer, The primary limitation on the use ofa an electric a 
computer appears to be the difficulty in modifying the circuits to reflect 
changes” in the agp of the structure. Each such change involves consider- 


‘computer, but at the end of the work described by the author all enh endl 
was gone. . There remained d only conviction that accurate and 


wer 


good and that those who operate the computer are © highly competent. 


1960, by R. H. MacNeal, (Proc. Paper 2578). 
22 Cons. Engr., Leeds, Hill and Jewett, Los Angeles, Calif. 


Ai _ _ The writer, who is responsible for the design of Blue Ridge Dam, was ni 
— 
~ 7 a 
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7 ‘computation of dynamic response of “ected with special application to male 
highway bridges. The author discusses the validity of his ‘assumptions _ 
- referring to experiments with single axle load. 7 Experiments \ with two- -axle 
loads are not yet reported in the literature known to the writer. However, at 
_ the Royal Institute of Technology, Stockholm, Sweden, a number of field — 
have yielded r results that, in spite of imperfect conditions, could give indica- 
7 tion of the validity of the theory on bridge structures. = | 7 
_ One of the tests was made on a girder bridge with 15.8 m (51.7 ft) span. 7 
. The loading was a two-axle truck with a total weight of 10.2 metric tons and 
a speed of 20 km per hr to 70 km per hr (12.4 mph to 43.5 mph). The vege 
surface was smooth in one series of tests but had corrugations of 3/4 in. -by- . 
in. boards in other cases. The dynamic moments and reactions were ‘computed — 
by equations given by the author. The differential equations were a 


In many cases— _the- magnitude of the dynamic increment of the bending 
oment at — in tests, 10% to 50%, cormingente s with computed values, 


mo 
an 
that 


« 


wise supports the assumption of the ‘shape of the ‘dynamic deflection curve > 
that, generally speaking, considers only the first mode of vibration. The rec- 

ote of the reactions at Supports compared with c computations indicate ee 
influence of higher frequencies, but the magnitude of the dynamic increm 

- can still be predicted by computations within a narrow a 

Opposite is the case in which “resonant” Speeds are considered. As 

could be expected, ‘the damping of the system can not longer be neglected and 
reduces considerably the amplitude of vibrations. Tests with “smooth” 
surfaces too give pooer a affinity to the theory | owing to deviations from —_— - 
a of the top surface on the bridge and the approach. a 

_ Acknowledgement. — The research program at the Institution of 

Engineering an and Bridge Building at the Royal Institute of Technology has been 
under the supervision of G. Wastlund. A. Hillerborg has been responsible for _ 

the test program and the writer has been responsible for the numerical com-— 

putations. The computer program was prepared by G. Hjerten. 


a October 1960, by Robert K. Wen (Proc. Paper 2624). | 
13 Civil Engr., Institution of Structural Engrg. and Bridge Building, Royal Inst. of | 
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‘Discussion by A. Assur 


loads, prepared by the author has been needed for a considerable time. The 
recent advances in the of ideally plastic plates have 
o ‘It may be argued that the stress distribution leading to the full plastic bend- 
‘ing: moment Mo = foh2/4 does not correspond to the true stress picture in an 
ice” ‘sheet. For fresh ice and quick loading this argument may | hold; however, 


& 

More” serious in its consequences is s in his Eq. 
‘that the hydrostatic reaction is confined within the hinge circle. There is no a 
reason why this should be so. Moreover, large scale loading tests by the 
laboratory with which the writer is associated have shown that a considerable 
deflection exists at the circumferential hinge ‘until collapse « occurs. After this 
+ ice rises at this point and only then is the hydrostatic reaction essentially 
confined to the hinge circle. Actual breakthrough occurs only with additional 

loading or after additional time. Meyerhof, neglects shearing forces and a re- = 
es on the surrounding ice sheet, which presumably acts as an 


ys Meyerhof overestimates the hydrostatic reaction, therefore his resulting 
radial moments at b are too small. 
tt is typical for solutions of this type that some of the geometry cannot be | 
_ predicted but has to be assumed. Meyerhof d does not predict the hinge radius b 
but, incorrectly, assumes it equal to the radius of the deflection dish in the 
_ etaatic case (b = 3.9 L). This is not supported by tests. The actual circum-_ 
_ ferential crack should and does occur well within the deflection radius, a fact 
_ Meyerhof discusses as proper but neglects in his further derivations. ee 
i hof is forced to stipulate b = 3.9 L because he assumed the full hydrostatic © 
_ reaction within the hinge circle to begin within Eq. 47, This is an inconsistency — 
’ that should be corrected in further work on the subject. Ver tna? 
There is good reason to assume that the circumferential crack should ae: 
at the place where the elastic theory would predict the maximum Pia 
oem. Tt can | be » shown that the following c condition holds for this location: 


a. October 1960, by G. G. Meyerhof (Proc. Paper 2627), 
Chf., Applied Research USA Cold Regions Research and Engrg. Lab., 
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numerical evaluation it was found that the result, can be 


: _ _ The definition of L by Eq. 98 is subject to some criticism. This expression _ 
is: based on the assumption of Hooke’s law with a a triangular stress distribution 7 
in the vertical profile. The plastic analysis, however, is based on a rectangular - 


t istributi mi t mo 
s ress distribution. Assumin ng furthe hermore w= of 


the “collapse of ‘on an elastic foundation, based o on the concept of 
yield lines. Johansson gives a rather derivation, containing errors. 
The writer, using the same approach, derived the following condition for i 


oe 


20 « “Forsok med armerade betongplatter pa under lag,” by A. 
(1947, 
tion), 


oe 


old quite well. Young’s modulus. J 
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DISCUSSION 


ave 


This may be compared with the corresponding approximate /~ 6) of Me er- 


— On- the basis c of Johansson’ s theory one can also predict the location of the 
circumferential er crack ick as bh 


in 


eye 


“This agrees 1 fairly well on | Eq. 98 i in a given range. For ‘small eoding 
areas the circumferential crack should approach the center rapidly, as is 
_ The “stiffness modulus” to be weed in the computation of lo 
Eq. 99 is very small. A value of 2300 sq tons, for example, was found for _ 
A further theoretical development along this line shows that the deflection — 
at the circumferential crack (before collapse) should be 3/4 of the maximum 
deflection at the center, while Meyerhof assumes | zero. Actual tests show 
’ ‘considerable deflection ‘at this” point although not qu quite 3/4 of the maximum. 
A number of tests also have been made by our laboratory to study the | 
punching of ice under concentrated loading. Meyerhof suggests use of Eq. 36 4 
7 for this phenomenon. He overestimates, however, the ultimate load, because 
Eq. 36 completely neglects the horizontal tensile stresses that are created 
= to bending. We prefer ” add the corresponding stress vectors at — 


Meyerhof quotes approximate equations for the elastic cases. 
primarily introduced by Westergaard, are no longer necessary, because very 7 


detailed | tables of the Kelvin functions vacated the exact solutions are — 

Throughout his paper uses many other approximations even 
though more precise be derived. This can easily produce 


-“ edge of a semi- is. 

in term of the unknown b. On the other hand 55 gives the 
collapse load for an infinite ice sheet. The ratio of Eq. 55 to Eq. 92 is then 
a convenient number : showing how much smaller the load has to be if os a 7 : 


a | “Tables of Kelvin functions and their derivatives, USA, Snow, Ice and Permafrost | 


Research Establishment, ” by D. E. Nevel, Technical Report 67, 
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ag Fig. Poa shows the poe for various a. For small e/b the curves and their 
first derivatives appear reasonable. The lines should not curve upwards, how-— 
ever, but downwards exponentially approaching unity. The ratio exceeds unity 
by x more than 50% in the region for which Eq. 92 was declared vane. — 


be: Meyerhof limits his analysis: to the discussion of ‘single eleeion loads or 

strip loads. For” practical purposes it would have been 1 very useful to also’ 
discuss the case of two circular loads at various distances (idealized aircraft 


load) or even four circular loads (wheeled vehicles), 
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a ‘PHYSICAL METALLURGY AND MECHANICAL PROPERTIES 
_ OF MATERIALS: BRIHLE FRACTURE® — 


i n by H. L. Sujata and Hans Conrad 


L. SUJATA, 34 M. ASCE, AND HANS CONRAD, 4 —In the section Atomic: 


"Mechanisms of Cleavage Fracture, the author has emphasized the role dis- 
_ locations play in yielding and fracture of these materials. While the author _ 
has presented the main features of the current thought on this subject, ‘the’ - 
writers would like | to supplement his treatment by elaborating on the relation- _ 
‘ship between dislocation behavior and the mac roscopic stress-strain diagram. 7 7 
Fig. 8 depicts, schematically, the sequence of steps in terms of dislocations © 
“that lead to fracture (ofa typical semi-brittle, dense, fully annealed crystalline 
a solid. (The term semi-brittle is used rather than brittle to distinguish crys-_ 
talline materials that show some plastic deformation before fracture , although 
it may be extremely small, from fully brittle materials such as glass.) Shown 
are the approximate stress — ranges in which certain dislocation phenomena ~ 
have been observed to occur. The phenomena listed in this figure have _ 
been observed completely by any One investigator on a single material, but 
rather represent a composite of information that has been obtained by a num- 
_ber of researchers on many semi-brittle crystalline materials. as 
Fig. it is ‘seen t that on application of a stress, the movement of 


“new” ‘dislocations first occurs (Tm). These dislocations are so defined 

7 cause they are not the intrinsic dislocations i inthe annealed material, but have 
been generated by scratching or by sprinkling a harder material (such as - 
‘silicon carbide), on the surface of the specimen.” 6,7,8 These new dislocations 

_ Move more ) coamey't than the intrinsic ones that are generally locked by impurity 


As the load is further mania multiplication occurs and additional dis- 
; locations are formed | (7g) giving slip bands. These dislocations pile up at 
some obstacle resulting in a a stress concentration tl that is proportional to 


1960, by B. L. Averbach (Proc. Paper 2686), 


3 Senior Engr., , Atomics Internatl., A Div. of North Amer. Aviation, Inc., Canoga — 


4 Senior Tech. ‘Specialist, ‘Atomics Internatl., A Div. of North Amer. Aviation, Inc. 


“Dislocation ‘Theory Applied to ‘Structural Design Problems in Ceramics,” by H. 


: Conrad and H. L. Sujata, paper to be published by the Materials Advisory Bd., Natl Z 
a “Academy of Sciences, “Symposium, Design with Brittle Materials,” December, 1960, 


aaa “Dislocation Velocities, Dislocation Densities, and Plastic Flow in Lithium Fluo- _ 
ride Crystals, ” 


by a © Ww. Gd Johnston - - Journal Applied Physics Vol. 30, Pp. 
+129, (1959). 


7 “Environmental Effects « on the Mechanical Properties of Ionic Solids with Particu-_ 


lar Relerenees to the Joffe Effect, ” by R. J. Stokes, T. L. Johnston, C. H. ais ‘Trans-  § 
AIME 218, p. 655 (1960), M7n 
— D. , F. ‘Stein and J. R. Low - J. App. Phys 


Vol. 31, 362 (1960). oe 
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FIG. 9. OF GRAIN SIZE ON THE LOWER 
-YIELD STRESS OF ELECTROLYTIC IRON. 
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7 number r of dislocations i in the pile up. up. On further increase in n stress, the stress ; 4 
os level at the head of the pile-up becomes | sufficiently large (as a result of the : 
stress: concentration) to give yielding (7 y) and the nucleation of microcracks 
(Tn). Because the number of ‘dislocations in the pile-ups is theoretically _ 
proportional to the grain size, the effect of grain on the yield stress 
given 6 is — of this relationship is 


-&# Prenagybe vs some | circumstances, be temporarily halted before final fracture 
occurs” (7,). In the region Tq to 7T¢, the effect of grain size on the fracture 
stress is given by Eq. 8. Support! Ofor this relationship is found in Fig. 10 for 
mild steel. It is noted from this figure that at large grain sizes (such as small — aa 
values of d-1/2), the fracture stress is equal to the yield stress, whereas at : 
ie “small grain sizes the fracture stressis larger than the yield stress and obeys 
_ the Griffith-Orowan relationship, the « author’s Eq. 2. Some e indication that 
ceramics also these fracture relationships11_ shown in Fig. 11. 
It is evident from the preceding that dislocation movement and multiplic 7 
_tion play an important role in the behavior of semi- -brittle m materials. ‘Fore ex- 
ample, if the dislocations do not move or multiply, fracture can not take place. 

This information might be useful in selecting a working stress level in design 
problems. design methods for semi-brittle materials generally (in 1961) 
strength data obtained from modulus of rupture | tests. These data are 
_ treated in any number of different ways ranging from the simple application of 

a factor of safety to the modulus of rupture to the more exacting application — ; ; 
presented by W. -Weibulls. 12 Perhaps another way of obtaining a working 
stress is to relate it to dislocation movement and multiplication. A working om 7 
_ stress level chosen as the stress at which dislocations begin to move or multi-— fal 
ply should eliminate the possibility of fracture or failure. These stress levels” - 
have not been evaluated (as of 1961) in materials of engineering usefulness and 


9 H. Conrad, G. Shoeck - Acta Met. 8, 791 (1960). 
“Properties and Microstructure, by J. R. Low - ASM (1954) p. 
- 11 “Dependence of Mechanical Strength of Brittle Polycrystalline Specimens on_ Por- 
_ osity and Grain Size,” by F. P. Knudsen - J. Amer. Ceramics Soc. 42, 376 (1959). B7c 
12 “Statistical Theory of Strength of Materials,” by W. Weibul , Ing. _Vetenskups i} 
Akad. Hand., 151, Ceramics Abstract 1940, 
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PHYSICAL METALLURGY AND MECHANICAL PROPERTIES 
OF MATERIALS: FATIGUE OF STRUCTURAL MATERIALS2 


Discussion by S. K. Ghaswala 

GHASWALA, 4 F. sack. one inand 


the progress and | development of aluminium alloys in structural engineer- 


ing since about 1935, , the writer would like to deal with a few salient aspects: 7 


_ The search for basic knowledge on fatigue that is to discover the mechanism | 
responsible for the formation of cracks, began over a century ago. The pro- 
‘gress of this search is evident from the growing mass of symposia and litera-__ 
ture that has slowly accumulated in the field (including the field of fatigue be- | 

havior of aluminium alloys.15,16,1 17, 18, ,19,20, 21,22,23 However, in spite ofall 
recorded advances, no satisfactory fundamental theory of fatigue relating the 
phenomena with the established Natural Laws has been formulated.24 This 
ean was - strikingly revealed in the two disasters of the British planes about 1954, 
Comet aircraft ALYP met with accident : at Elba, and Comet 
G-ALYY crashed near Naples. A Court of Inquiry was appointed to go into the © 
details of this accident. One of the three Assessors was WwW. Duncan, who 

carefully reviewed the causes of failure and gave ana admirable survey of he 
‘aps in our extant knowledge onthe subject.2500 


a — to Duncan one of the unexplained phenomena of fatigue in metals 
s size effect, | observed in light of the fact that the fatigue strengths of geo- 


| metrical similar specimens of the same metal are not proportional to their 


_ 15 “Mecanique Industrielle,” by J. V. Pongelet, 2nd Ed., 1839. pas 

16 “Resistance des Materiaux,” by A. Morin, 1853. 
= “Internal Stresses and Fatigue in Metals,” 

«18 “Fatigue in / Aircraft "Structures, by A. Academic N.Y 


December 1960, by Horace J. Grover (Proc. Paper 2688), | 
1 


1956. 


19 “Metal Fatigue,” edited by J. A. Pope, Ch. and Hall, London, 1959. per ce 
_ 20 “Fatigue of Metals and Structures,” edited by H. J. _ Grover, S. A. Gordon, and L. 
Jackson, Thames and Hudson, London, 1956. 
_ 21 “Bibliography on the Fatigue of Aluminium Alloys,” by D 
q over 480 references), Aluminium Dev. Assoc., London Tech. Memo 147.” 
22 “Fatigue, a Bibliography,” by J. L. Erickson, Light Metal Age, November, 1945, 
28 3 “Fatiuge of Aluminium, ty Templin, (contains nearly 500 references), 
‘Proceedings, ASTM, Vol. 54, 1954, p. 641. 
Some Aspects of the Plastic Design of Aluminium Structures,” by S. K. Ghaswala, 
Assoc. for Bridge and St. Engrg. Pub., Zurich. Vol. 16, 1956, p. 231, = 
_ 25 “The Comet and Design against Failure,” ” by W. J. Duncan, Engineering, February _ 
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cross sectional areas—the safe stress actually being re reduced with increase 


_ The phenomena of stress concentration, or as it could be rightly referred 
to “stress intensification,” as suggested by Duncan, playsa significant part in 
fatigue and plasticity of metals. Inthe case of aluminium alloys and some other _ 
metals which do not have a well defined yield point the plastic flow when 
stressed beyond the elastic limit permits the relief of high stress concentra- 
at the highest However, in 1adjacent regions the stress increases. 


nature of loading and intensity, and type of material used, it c can be said that — 
a slight amount of plastic deformation is capable of ‘relieving appreciable 

stresses. . The result of plastic deformation in a member having holes and 
cutouts or any other form of stress raisers, influences its behavior under re- 
peated loading and may increase its fatigue life considerably. The overall de- 
sign of complex structures for fluctuating loads therefore assumes complica- 
tions unlike steel where a definite fatigue limit exists. 

_ Stated in engineering parlance as understood by civile engineers, aluminium — 
alloys display differing fatigue properties | with differing alloys and strength. 
In high-strength alloys there is not much difference i in the basic fatigue proper- 
ties of alloys of the type of Duralumin and the Al- Zn-Mg- Cu alloys. With 


sults. in n appreciable difference in 1 the fatigue strength, although the alloys. may 
_ have virtually the same tensile strength. Compared with the heat-treatable 
_ alloys o of aluminium, the Al- ‘Mg group shows a strikingly sharp k knee in the 
- fatigue ¢ curve and a definite limit in certain cases. In the case of low- strength 
_ alloys, including commercial purity alt aluminium, the fatigue curves tend to 
a typical test using aluminium alloys having tensile strengths ranging 
from 5.5 to39 tons per sqin. it was found that the fatigue ratio (ratio of fatigue _ 
_ Strength to ultimate tensile strength) ranged from 0.26 to as high as 0.59. AD a 
i of these results indicated that for non-heat-treatable alloys the fatigue 
- ratio lay in the higher range from 0.45 to 0.59, and could be compared with 
corresponding ratios for mild steel. On the other hand the fatigue ratio for - 
_ heat- treated materials was very much lower, indicative of doubtful fatigue 7 1 
limit. These figures were based on an endurance of 50 million cycles, which 
it must be realized is a comparatively very high figure, because in practice 
_ this is equivalent to one -eycle per minute for 95 yr, or 20 cycles per hr for 
285 yr. Normally 50 loadings (alternating) per working sateen build up to 
300,000 loadings in approximately 20 yr. 
Asa broad generalization it can be saidthat as the fatigue ratio drops with | 
_ increase in tensile strength, ‘it becomes obvious that if a structure is designed | 
on a purely static strength basis its fatigue resistance will lessen as the static a 
strength of the material increases, ie - Furthermore for static structures in 7 
Boek which carry dynamic loads, such as bridges, the dynamic- static 
load ratio would invariably be greater than for a similiar structure in steel, 
because of the lighter weight of the structure itself. Fatigue problems tnere- q 
fore become more accute. The case of high-strength materials, especially _ 


those having a very high ratio. of yield stress (or proof stress) to ultimate — 
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"stress, requires special of | stress concentration in such 
—_ applies to materials other than aluminium alloys— the only | 
real limit to the stress which may be in the material is the limit of propor- 

a tionality or the. ae. stress. In other words it means that the stressing in in 


‘mean stress in the « case of a a . material having a a high yield stress than in one 
having a low yield stress. As such, apart from fatigue characteristics, = 


considering the superimposition of a reversing stress+ z on a steady stress 
condition Fs. If Fo is” the static failing stress and F, is the fully reversed 
stress which would cause the member to failat a given number of cycles, then 
it can be assumed that in the condition representing failure eo 
= 
It is therefore essential to k to —_ thooen sum — the stress ratios below unity while a 
in no case should Fg + Fy exceed the allowable design stress. 
= The civil engineer has much to learn from the aeronautical designer since ; - 
‘an aircraft forms one of the most useful and time-tried indices of gaging and 
solving fatigue problems. This is evident from the fact that during some 40 yr 
of failure investigation by C. W. George, 27 over 75% of the specimens examined . 
failed by fatigue, giving a crystal clear indication of the importance of this © 
_ Thus, the methods of providing fatigue safety in aluminium includes the 
Ss adoption of rational engineering principles by avoiding | stress raisers, and 
gaging, as best as possible, the actual alternating loading which would come 
on the structure during its lifetime. No amount of theory can supplement — 
simple fact that absence of a small stress raiser is far more important than — 
mere strength, and that a simple butt weld is superior in fatigue to the same 7 


weld heavily reinforced by as late (or a rib. 


26 “The Strength of Aluminium,” Aluminium Co. of Cani Canada | Ltd., 1960, 


dournal Birm rm Met. ‘Soc. Vol, 2 27, 1947, Pp. 308. 
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ELASTO-P PLASTIC ANALYSIS BY BY) NUMERICAL PR PROC 


Discussion by Le-Wu 

A. M. ASCE. —The author has developed some ‘useful methods 

for - determining the deformations of partially plastic beams and the ultimate © 

strength of eccentrically loaded columns of aluminum alloy. In a recent invest- 

igation9 on inelastic instability of frameworks, the writer has adopted the 

numerical integration technique devised by von Karman10 (that is also the 

basis of the scheme shown in Table 1 of the paper) to study similar problems. __ 

To determine the deformation of a simply supported beam ‘subjected to 

elasto- plastic loading, it is “only necessary to perform a series of direct 


integrations, segment by segment, using the predetermined moment- curvature 
(M -¢) relationship of the cross section. For the particular type of stress- 
aon egy assumed by the author, the M- ¢ relationship may be obtained 


22 canbe used to solve for or M/M, to various assumed values 

. of 6/dy. y- Results of computation for three B values are shown in Fig. 16. Using 

these hhon- ‘dimensional M @ curves the deflected shape of the beam of 

: ample 1 can be . obtained by integration, _ starting from the center of the span 
j in which the slope is known to be zero. The integration procedure is similar _ 

to that illustrated in Fig. 13. In Fig. 17 the center deflections thus determined | 
= plotted as the dot-dashed line against the integration interval a chosen in 7 

the computations. (This is to study the influence of A on the accuracy of the 


tha following may be observed: nial 
The integration “interval has considerable influence on the results 
7 a; 2. ‘Newmark’ s procedure yields an upper r bound to the true deflection. ‘The 


or for this is th that | the formulas used to estimate t the equivalent concen- | 


» 8 Research Assoc., Fritz Engrg., Lab., Lehigh Univ., Bethlehem, Pa. a: = 
7 9 “Stability of Elastic and Partially Plastic Frames, ™ by Le-Wu Lu, Thesis presented © 
7 to Lehigh Univ. at Bethlehem, Pa. in 1960 in partial fulfillment of the requirements f 
10 “Untersuchungen Uber Knickfestigkeit, by Th. ¥. Kirmdn, Forschungsarbeiten, 
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trated loads in such procedure — seit give higher values for cases in 


which the M/E I diagram has sharp variations. 


-- 3. The method of direct integtation always leads to a lower bound to the 
actual deflection. This is peeneny due to the under-es estimation of of the the average 


Thus by averaging the results obtained from these two procedures, a better 
Ee of the deflected shape | can be achieved. The solid line of Fig. 17° 
shows the average center deflection of the beam under consideration. It may — 
a be noticed that even with A = 6 in. the center deflection can be determined by 
a averaging the upper and lower bounds of the solutions t to within 6% of accuracy. 
_ Referring to the problem o of instability of eccentrically loaded columns, — 
the author’s approach seems to be similar to that developed by von Karman 
discussed ‘Timoshenko. 5 In such an approach a set of column de deflection 
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16 LASTIC MOMEN T-CURVATURE RELATIONSHIPS: 


curves is is: first constructed for a specified value of the axial force P. Each 


4 

elo ongation of the two extreme fibers of the height section. The critical 
length of of a column can then be determined by the graphical method indicated _ 
‘Fig. . 15. Morris Ojalvoll M. ASCE has developed a more systematic 
n method for determining the critical length of steel columns, in which — 
column deflection curve | is constructed for specific value of the initial 
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> DISCUSSION 


“Fig. 6 will be described ‘to ron ‘the manner in which Ojalvo’ s s method can “— 
applied to columns of aluminum alloy. The particular problem selected is as 
follows. 


Given: (1) Aa Axial force P=0.4 fe A, in which A is the cross sectional area 
of a rectangular column; (2) ) Eccentricity € = 1.0 R2/c, This corresponds to 
an end moment of M/My = 0.4. (3) The stress- _strain diagram of the material 

_ Required: The maximum column length that is in n stable equilibrium | with 


the applied axial force and the end moment - 
‘The non- dimensional - M - A diagram of Fig. 12 may be interpreted as the 


by half (because A /ey = 20/dy. ‘The curve corresponding to ) F/fy = =0.4_ 
can then be used to construct the required column deflection curves by the | 
numerical integration procedure (Table 1), starting from an initial point in 

_ which the slope is equal to §,. Severalsuch curves are shown in Fig. 18. The © 
_ information from these coleman deflection curves is then organized as shown © 
in Fig. 19. In the upper portion the moment at a point on the column curve is 
plotted against the value of the slope 6 at that point. In the lower section the — 
distance of the point on the column curve from the mid-height (non- dimen- 
sionalized as L/2R) is plotted against 4. A horizontal line in the upper portion — q 
is drawn at the specified moment value. Intersections with the M - 4 curves 
are carried down tothe corresponding curves inthe lower section. This deter- 
mines the equilibrium lengths for the selected column configurations. The 
longest length found gives the solution. For the present problem the maximum > 
slenderness ratio is equal to 48.2. It may be seen from Fig. 6 that this value of © 


L/R is in s in agreement with the author’s result. 
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oy WIND I INDUCED VIBRATIONS IN ANTENNA MEMBERS2 


Discussion by Glenn B. W Woodruff 


B. WOODRUFF, 3 ASCE.—Reliable data for the von Karman coeffi- 
- cient ar are needed not only for antenna members but also for larger construc- 
_ tion such as stacks and refinery vessels. It is to be hoped that the author con- a 
_ tinue his 3 experiments and furnish such | data through the entire range up to 
“a Possibly, the author’s presentation would be clarified bya better definition : 
of the two cases, the first in the Reynolds sub-critical range, that of forced 
vibrations with the rate of vortex sheading governed by | the Strouhal number 
and, the second, that of self induced vibrations with the rate of vortex sheading © 
governed by the natural frequency of the vibrating member. These two cases 
are well illustrated by. Peter Price , (23), Fig. 7. It may be noted that in n the 
"second case, there is no “critical” 
_ The author fails to explain the mathematical processes by which = derived 
his values of c.. In their analysis of the failure of the Michigan stacks, with 
an R of 7,000,000, Dickey and Woodruff, (24) assuming a magnification factor of 
1/6, checked their assumption of Ck = 0. 66. If a factor of 7/5 were assumed, 
im & "2 would become 0.21. This is in close agreement with Ozker and Smith (25), 


7 who arrive at the value of large stacks of 0.20. These two investigations indi- 
cate that C, is constant above R=500,000, 
A factor that has received little investigation is that the wind induced ampli- 
‘tudes, rather than reaching a steady” - state, show 2 a tendency towards beats. 
This is indicated by the authors Figs. 2(b) and 4(b), Penzien’ s Fig. 4, Fig. 3. 
of | Ozker and Smith (25) and numerous other oscillograms. F. - ‘Farquarhson 
_ F. ASCE, (26) found that the frequency of wind induced vibrations was consist- 
- less than the natural frequency, the difference increasing with the damp- 
— ing. _As against Farquarhsons’ observations of lesser frequency under wind, 
Dockstader, Swiger, and Ireland (27) report, for the Morro ‘Bay unlined stacks, 
__ observations confirming the tendency towards beats but a a wind- induced fre- — 
quency of 73 cycles per minute, higher than the natural frequency of 67 cycles 
A further - investigation « of this matter may y show that the response is not as 
7 ‘random: as is indicated cialis author and that, with appreciable damping, the a 
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ied in phase are not sufficient to ‘reach the magnifica- 
check of the numerical value o of R in Appendix I. : 


The writer suggest 
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23 “Suppression of the Fluid Induced Vibrations of Circular Cylinders,” ” by 

Peter Price , Proceedings, ASCE, Vol. 82, October, 1956. 

by Walter | L. Dickey and and Glenn B. -Wood- 


O pity Vibrations of Steel Stacks,” 

ruff Transactions, ASCE, Vol. 121, 1956, p. 1054. ; 

29 “Fa tors Influencing the Dynamic Behavior of Tall Stacks Under the Action. 

att of Wind,” by Ozker and Smith, Transactions, ASME, Vol. 78, No, 6, August, 


1956, p. 1381. 


“Which Influenced the Design of ‘the New Tacoma 
’ by F. B. Farquarhson, Univ. of Experiment 
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Williams, October 1960, Proc. Paper 2619. 

«Pp. 3 36 and 37. 37. In Figs. 2 ar and 3 3, change cHtoC 

41 and 42. In the captions for Figs. 4 to x 


43. In the fourth line below Eq. 42, change Figs. to Eqs. 


59. In the line Bq. 4b delete. supers¢ cript 
Pp. 62. In line 1 0 1 of f paragraph 5 delete the w wor rd were 7 


?p. 64. In the line below Eq. 15 change Hin( ) to Hn (1 ) ). 


( 
In the tl third line below Eq. 15 change Tol )(D) to Jol da 
In the fourth line below Eq. 15, change" Hj )t 
p. 66. In Eq. 21, change the term m (1/ 73) w to | (1/r 


mw 
74 17. In the c: Figs. and 8, “change V to 


76. In Eq. 73, =" v vt in three 


3.73 Mp = 0. 


and for P9 as follows: 


0.928 Py 
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i 
ange Jj() to 
= 
= -47Mp / | loge +k| = 7.46 Mp = 0.377 Pu 
80. Rearrange the expression for P3 as follows: 


p Pp. . 88. Eq. 104 should be written 2 as 


“Dynamic Response of Beams Traversed by Two-Axle Loads” by Robert K. : 


Pe 93. _ In line 2, delete the parentheses a around the words “an infinite num- 


ber r of degrees.” 
p. 95. In Eq. 3, change th the term mg to ) mp. In Eq. 7, /, change the term f (t) 
fo). 


96. In Eq. 9, change Z to z in two places. 


101. In line 4 of | the ‘seventh full paragraph, t/L : =0, 55 to to t/L 


October 1960, Proc. Paper 2624. 


In line 2 of the ‘fifth full peengrem, change the reference number 7 7 


. 105. In} line 2 of the last delete the word equal. 

. In line 5 of paragraph 4, change the reference neaber 7 to 8. 

of paragraph 6, change the reference number 12 to 13, ae 

107. » the expression for Ag, in the first line, change the order of ap- 


no of the final parentheses and the “subscript j. In the expression for 


insert the term J" just after = and before the subscript j. =~ 
108. In the the ‘first line (of Eq. 41, insert a bar over the term ey. 


line of Eq. 41, insert a immediately after ‘the 
term 1 Ee In the first line of Eq. 42, insert a bar over the term 5? . 
109. ih the after Eq. 43, insert a bar over the term 


In the line above ™ 45, insert t the t term, n after the word in and before the: 


“111. Inth of vj . insert the term immediately 


In Eq. to 


In line Ea. 4=2.2Ltor 


and for Pq as follows: | 
= 16M 1 - = 19.75 My = P 
lt 
q 
: 104. 
ex 
a 
— 
I 3 — 
— 
ae —?p. 110. In the definition of €j insert theterm ¢€ immediately after ; and be ries 
— 
960, 
— 
2.2 


word Infinite to the solid curve. | 


p. ‘118, In Fig. 3a, add “Infinite:” 2 
122. In line 15, change v2 2a to 2a 
lin 


123. In Eq. 25, ‘change a V2 tos 


p. 124, “In the first line after Eq. 26, change 2 VE Lt to 2 Vet. 


. 131. In Fig. 8, 8, add “in, to the abel “Values fp.” 


line 3 bottom, p. 135, line 5 of Eq. 43, p. 137, ‘first. 


after Eq. - = 


op. 136. In the numerator of Eq. 50c, insert bracket | between n the terms P 


and and another bracket just after the term b2, 


ip. 198, Int the fourth line after Eq. 62, insert < betw en the terms a/b and 


noi | by Morris Ojalvo and Ahmet Cakmak, October 1960, Proc. eon 
2632, of “ ‘Behavior | of Buckled Rectangular Plate” ’ by Manuel Stein, April ened 


Pro Paper 2445, 
— 


-p.170. In Fig. 16, Ref. 18 to Ref. 20. 
‘teaciain by R. J. P. Garden, October 1960, Proc. Paper 2632, of “Commen-_ : 
ta tary on Plastic Design in Proc. Paper 2454. 
theses should be added immediately b before and after the minus sign preced 
“ing In Eq. bracket should be inserted immediately prior tothe 
‘mins ane caged 6'' 34. Just after the same minus sign, a parentheses 
“should be added. The “bracket t following the t the term 6" "348 should be changes 


a parentheses. 
179. In Hine 1e 5, chi vature” to “curved.” 
oP. , 179, In line 5, ¢ ange “curva ure” to “curve 
Fixity Effect on Vibration and Instability” by Avid Burgreen, December 


e. 19. The equation presently denoted as Eq. 11 is — the first half of 
Eq. 11, the last part of which was erroneously labeled Eq. 12. The 


we 
In the third line above Eq. 1, change @ nt 67 In Eq. 
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